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We study the scissors mode (orbital Ml excitations) in small Na clusters, triaxial 
metal clusters Nai2 and Naig and the close-to-spherical Nag"*", all described in DFT 
with detailed ionic background. The scissors modes built on spin-saturated ground 
and spin-polarized isomeric states are analyzed in virtue of both macroscopic collec- 
tive and microscopic shell-model treatments. It is shown that the mutual destruction 
of Coulomb and the exchange-correlation parts of the residual interaction makes the 
collective shift small and the net effect can depend on details of the actual excited 
state. The crosstalk with dipole and spin-dipole modes is studied in detail. In par- 
ticular, a strong crosstalk with spin-dipole negative-parity mode is found in the case 
of spin-polarized states. Triaxiality and ionic structure considerably complicate the 
scissors response, mainly at expense of stronger fragmentation of the strength. Nev- 
ertheless, even in these complicated cases the scissors mode is mainly determined by 
the global deformation. The detailed ionic structure destroys the spherical symme- 
try and can cause finite Ml response (transverse optical mode) even in clusters with 
zero global deformation. But its strength turns out to be much smaller than for the 
genuine scissors modes in deformed systems. 
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I. INTRODUCTION 

The scissors mode in nuclei is a small-angle collective rotation of a deformed proton cloud 
against the complementing deformed neutron cloud. It thus belongs to the basic magnetic 
orbital dipole excitations. A rotation-like character of the mode assumes that it exists only 
in the systems with broken spherical symmetry. Early theoretical predictions are found in 
Refs. 1^ and its first experimental observation was reported in Ref. 0. For a review 
see Ref. Since then, the scissors mode has also been discussed in other deformed 

systems, for example in quantum dots P, Bose-Einstein condensates 0, and last not least 
in metal clusters [0, ^. In the case of deformed metal clusters, the mode is viewed as a 
rotation-like collective oscillations of the electron cloud against the ionic background. There 
exist already a few theoretical considerations on that topic, early estimates from a collective 
perspective [^, RPA calculations with a phenomenological Woods-Saxon potential P], and 
fully detailed RPA calculations using density-functional theory and jellium for the ionic 
background 0. The relation of scissors modes to orbital magnetism had been discussed in 
Refs. [0, H, H, |ll| . Experimental access to scissors modes in clusters is yet an unresolved 



problem ||T2| (see also discussion in Ref. 0). Nevertheless, or just because of that, further 
theoretical investigations are worthwhile in order to elucidate the nature of these modes and 
their key features. 

It is the aim of this paper to present a theoretical investigation of scissors mode in Na 
clusters using time-dependent density functional theory for the electrons and fully detailed 
ionic background. We will consider in particular the triaxial clusters Nai2 and Naie. These 
two test cases provide several interesting features. As triaxial clusters, they have deforma- 
tion in all three spatial directions and thus should display three distinct scissors modes. 
On the other hand, the single-particle excitation spectra are most diffuse in these soft sys- 
tems which may give rise to strong fragmentation. Both clusters have spin-polarized and 
axially symmetric isomers which give rise to crosstalk with the spin-dipole mode which is 
energetically very close to the scissors mode. (The dipole and spin-dipole excitation spectra 
in Nai2 had been studied extensively in Ref. [11^.) The axial isomers allow also a direct 
comparison between scissors modes in axial versus triaxial shapes. We will investigate the 
impact of ionic structure on the scissors modes in two ways: first, we compare the triaxially 
deformed test cases with the results of the smooth jellium model, and second, we consider 
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as a complementing example the nearly spherical Nag with full ionic structure. The test 
case of axially symmetric Nan'*' with jellium background is used to disentangle the various 
(counteracting) parts of the interaction, which compose the final prediction of the frequency 
of the mode. 



II. FRAMEWORK 

The electron cloud of the clusters is described by density-functional theory (DFT) at 
the level of the time-dependent local-density approximation (TDLDA) using actually the 
density functional of Ref. fl^. The coupling of ions to electrons is described by a local 
pseudo-potential which has been proven to yield correct ground state and optical excitation 



properties of simple sodium clusters |T^. The electronic wavefunctions are represented on an 
equidistant grid in three-dimensional coordinate space. The ground state solution is found 
by iterative relaxation. The time evolution is done by the time-splitting method. For details 
see the review |]16[ . 



To compute the spectral distributions, we employ spectral analysis after TDLDA prop- 



agation ||16| , p77[| . This technique which had been developed and much applied for studying 
dipole resonances 0, |l8l is equally well adapted to the present case of the scissors mode. 
The dedicated initialization of the scissors channel is achieved by rotating the ground-state 
electron cloud by a small angle out of the equilibrium. The subsequent dynamics is ana- 
lyzed with the orbital angular moment < Lj > which is the appropriate observable for the 
orbital Ml strength. Crosstalk between the principal axes of the cluster is very small. Thus 
rotation about the x-axis (using as a generator) explores the dynamics of L(t) = (L^) 
and similarly for y and z. We follow the time evolution for 200 fs. The Fourier transform of 
L{t) provides the spectral strength distribution of orbital Ml strength. We also study the 
crosstalk to other modes. This is done simply by taking a protocol of the dipole moments 
Di oc Tj as well as spin-dipoles Si oc rjO"o (where (Tq is the Pauli spin matrix along quantiza- 
tion direction) and Fourier transforming them. In that case, we consider the spectral power 
distribution which is the sum of squared absolute values of the respective Fourier transforms. 
For a more detailed discussion of the spin-dipole see Ref. |13|. 



The global deformation of the cluster ionic background is crucial for the scissors excita- 
tions 0, §]. It is characterized by the dimensionless quadrupole moments recoupled to the 
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total deformation (3 and triaxiality 7: 

/9= V^2o + 2/5|2 , 7 = atan 



72/^22 
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where iV/ is the number of ions (atoms). The definition of the (32m contains a simple classical 
averaging over the ionic positions. 

In order to check the impact of ionic structure, we also use, as an alternative to the 



detailed ionic background, the soft jellium model []T9| for the ionic density 



' l + exp{{r -R{Q))/a) ' ^ ' 

R{e) = Ro{l+ E S2mY2m{e)] (3) 
\ m=0,±2 / 

1/3 

where Rq = CrgNj , = 3.96 ao is the Wigner-Seitz radius, the coefficient C is adjusted 
to ensure volume conservation / drpi{r) = Nj, pjQ = 3/{AiTrs) is the bulk density. The 
diffuseness a of the jellium surface allows to achieve a good reproduction of the empirical 
optical properties |]19|, It can be justified by folding a steep jellium drop with a local 
ionic pseudo-potential [|^. We choose the diffuseness a = 1 ao. The generating deformation 
parameters ^20 and ^22 = <^2-2 are adjusted to give the intended deformations (3 and 7 
according to Eq. ([l|). The averaging in Eq. then becomes a weighted integration over 
the jellium density p/. 

III. BASIC ASPECTS OF SCISSORS MODE 

Before proceeding to the fully detailed TDLDA calculations, it is worthwhile to look at 
macroscopic features of the scissors mode as well as to present its simple microscopic view 
in terms of a deformed shell model. 



A. Collective picture of rotational vibrations 

In a collective picture, the scissors mode can be viewed as small-angle rotations of the 
spheroid of valence electrons against the spheroid of ions (see Fig. 1 , part a)). More 
precisely, the displacement field of this mode consists out of a rigid rotation of the electrons 
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against the positive ionic background plus a compensating quadrupole term which serves to 
tune vanishing velocity at the surface M (see Fig. 1 , part b)).: 



uir] 



n X f+(3{l + (3/3)-^V{yz 



(4) 



The distortion of the momentum Fermi sphere generates a restoring force leading to rota- 
tional oscillations which constitute the scissors mode 0]. 

The mode is characterized by low-energy oscillations with strong magnetic dipole transi- 
tions to the ground state. In axial clusters, the mode is specified by the states \A'^ = 1^ > 
where A is the projection of the orbital moment onto the symmetry axis z and tt is the space 
parity. The corresponding energy and magnetic strength are [0, §] 
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B{Ml) = 4(1+ I I 0)^6 
2 



where Np is the number of valence electrons (we use here natural units mp 



(6) 
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The value B{M1) stands for summed strength of the degenerated x- and y-branches. The 
z-branch has vanishing strength for symmetry reasons. It is worth noting that B{M1) does 
not depend on and so is the same for different metals. 

In triaxial clusters, the scissors mode splits into three branches with the frequencies [0 
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and the strengths 
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where vectors run over i = x,y,z components. The estimates show that both the frequencies 
and strengths are proportional to deformation f3, which reflects the fact that collective 
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scissors modes can exist only in deformed systems. This feature may serve as an additional 
indicator of quadrupole deformation in clusters. The scissors mode has been already observed 
in several systems with different physical nature (nuclei, Bose condensate, quantum dots). 
This means that the mode is a general dynamical phenomenon persistent for deformed two- 
component systems. It is interesting to note that there is another kind of universal orbital 
magnetic mode, the twist M2 mode, which has been recently predicted in clusters . The 
twist may exist in systems of any shape and it is the strongest multipole magnetic mode in 
medium and large spherical alkali metal clusters 



B. Ml transitions in the deformed shell model 



In order to see which single electron states are the most relevant in building the scissors 
mode, we check the structure of the transition in the framework of the Clemenger-Nilsson 
model |2^, The major shells Af of the spherical oscillator remain a useful sorting scheme 
for slightly deformed clusters. Taking into account that A/" is connected with the space parity 
of the levels inside the shell as vr = (—1)-^, one may immediately conclude that A'^ = 1+ 
scissors mode has to be mainly generated by low-energy AA/" = transitions inside the 
valence shell and high-energy AM = 2 transitions through two shells. 

More information can be obtained from the expression for the orbital Ml transition 
element which in axially deformed clusters has the form 

{^,\Ui\^h) = IK,,.Ja,,a,+i (9) 

$:<L«nLV'^(^+l)-A.(A.+ l) . 
nL 

Here, the wave function of a particle (z/ = p) or hole (z/ = h) deformed state '^lJ=^fn2A = 
J2nL (^nL^nL{r)Yi\{Q)xi/2u IS Written as a superposition of spherical (nL)-configurations in 
which n stands for the number of radial nodes; labels the projection of the principle 
shell number JV = rix + Uy + Uz into z-axis. Eq. shows that the scissors mode is 
generated by Ap = A/^ ± 1 transitions between wave function components belonging to one 
and the same spherical (?T,L)-configuration. In spherical systems (nLA)-states belonging to 
such configuration are degenerate while in deformed systems these states are energetically 
split and so may be connected by Ml transitions with non-zero energies. This is just the 
origin of the scissors mode. Obviously, the energy scale of the scissors mode is determined 
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by the deformation energy splitting and so is rather smalL This explains predominantly 
low-energy {AM =0) character of the scissors mode. Just the low-energy branch exhausts 
most of the scissors B{M1) strength. The high-energy (ZVV =2) branch of the mode carries 
much weaker strength since the particle states involved into (AA/'=2) transitions include 
only small fractions of (nL)-configurations from the valence shell. 

One may show that in the harmonic oscillator space the Ml transition matrix element 
is proportional to the matrix element of the quadrupole operator r'^Y2i. This indicates a 
coupling between the high-energy (AA/'=2) branch of scissors mode and electric Xfi = 21 
quadrupole plasmon in deformed clusters p. 

Note that in spherical clusters the Ml transitions connect degenerated states from one and 
the same (nL)-configuration. These transitions have zero energy and, therefore, the scissors 
mode in spherical jellium clusters is absent. The ionic structure destroys the spherical 
symmetry even at zero global deformation /? = and thus allows some L±i -response in 
such clusters. It remains to be seen how strong this effect is in actual calculations. 

Eq. (^) hints that the strength of the scissors modes increases with the orbital moments 
L which are involved. Besides that, low azimuthal quantum numbers A favor the mode. In 
heavy clusters high orbital moments are accessible and thus these clusters can exhibit strong 
orbital effects ||^. Already at A^^e ~ 300 the magnetic strength B{M1) can reach impressive 
value of 300-400 /i^. The RPA calculations confirm the trend for the frequency but 
reveal in light clusters considerable fluctuations around the trend (|^) for the strength |^, Q . 
These fluctuations are due to quantal shell effects. Finally, it is worth to mention that the 
scissors mode determines van Vleck paramagnetism and leads to strong anisotropy of orbital 
magnetic susceptibility R, 13, 0, . 



C. Preparatory example 

As a first example, we consider the simple case of a strictly axial cluster. To this end we 
chose Naii^ with a soft jellium background (generating deformation ^20 = 0.38, 7 = 0). The 
symmetry allows only degenerate x- and y-scissors modes. The associated spherical valence 
shell (2s,ld) consists of 4 levels: Un,K = 220, 211, 202 and 200. In Naii + , the Fermi level 
220 is fully occupied and the other three levels remain empty. The only possible scissors 
transition is 220 -^211. The low density of particle- hole [Iph) excitations reduces fragmen- 
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tation of the Ml strength and thus one has only one pronounced scissors peak, the ideal 
scenario for a first exploration. The upper panel of Fig. 2 shows the spectral distribution 
of Ml strength. The full line is the result from the linearized TDLDA calculations p, |21| 
with full interaction. We use this clean test case to identify the various contributions of the 
residual interaction. The ground state is always computed with the full energy functional. 
But for the excitation spectrum, we consider also the terms of the residual interaction sepa- 
rately. One finds in the upper panel the spectrum for the cases of no residual interaction at 
all (unperturbed Iph response) and of pure Coulomb. The Coulomb interaction is strongly 
repulsive (like for the dipole mode) and results in a strong blue shift of the unperturbed 
strength. The blue shift for the total residual interaction is much smaller. This is due to 
the strongly attractive residual interaction from the exchange-correlation part of the func- 
tional. In fact, the exchange- correlation interaction alone produces so much attraction that 
the scissors mode acquires an imaginary frequency. The mode becomes unstable and is thus 
missing in the figure. 

The lower panel of Fig. 2 shows as a complementing information the radial profile of 
the residual interaction. The density dependent operator Q2i{r) of multipolarity A/i = 21 is 
in fact the (radially averaged) response potential —{d^E/dp'^)6p where 6p is the transition 
density and E is the energy functional. The plot shows the full residual interaction as well as 
its two separate contributions. Repulsive parts are negative in this representation. One sees 
nicely that Coulomb is strongly repulsive, exchange-correlation is attractive and the total 
effect is a moderate repulsion. After all we see that the net residual interaction is composed 
from two counteracting pieces. This inhibits any general statement about the expected shift. 
Blue shift as well as red shift can emerge depending on the actual transition density. But 
one can predict that the total residual interaction will always be rather small (with possible 
exceptions in special cases). 

IV. RESULTS AND DISCUSSION 

A. Nai2: global deformation versus ionic structure 

Figure 3 shows for Nai2 the scissors strength, i.e. the Fourier transformed signal Li{t), 
as computed from real time TDLDA (section II). The spin-saturated triaxial ground state 
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and spin-polarized axial isomer are considered with both ionic and jellium backgrounds. 
Triaxiality mixes the states with different A but still conserves the space parity. The ionic 
background makes the states parity mixed as well. 

The spectra for the triaxial ground state (left lower panel) show the modes Lx-, Ly, and Lz- 
The z-mode has one dominant peak at the low-energy end. Two other modes are strongly 
fragmented. A quick glance at the lower insert explains this. The density of Iph states is 
very large in this neutral and triaxial cluster and thus spectral fragmentation is very likely. 
The lower end of the spectrum coincides with the lowest Iph excitation which hints that 
the net effect of the residual interaction is rather weak. A simplified axial treatment can 
roughly explain the origin of three main peaks in x- and y-responses: in Nai2 the Fermi level 
is one-half occupied level ATn^A =211 and within the valence (2s- Id)- shell there should be 
three main scissors transitions: 220 —>■ 211, 211 202, and 211 —>■ 200. 

Note also the Ml strength around 3 eV. This is the high-energy branch of the scissors 
mode. It emerges from the scissors transitions through two quantum shells {AJV = 2) as 



mentioned in subsection |III B| . This branch exhibits a coupling to the A/x = 21 component 
of the quadrupole plasmon [H. The energy range coincides perfectly with the position of the 



quadrupole plasmon around 3 eV (it is placed approximately at yQ/5 times the frequency of 
the Mie plasmon which is here around 2.8 eV). The strength is strongly fragmented because 
the weakly bound triaxial Nai2 has a high density of Iph states in this energy region (actually 
not shown in the figure). 

The upper left panel of Fig. 3 shows results for a comparable jellium background. For this 
aim, we have computed the global deformation of the ionic configuration. This yields (3 = 
0.55 and 7 = 17°. The jellium background is tuned to have precisely the same deformation. 
The Iph states start now at much lower energy since the energy gap at the Fermi surface 
is lower than with an explicit ionic configuration. This has two reasons: the jellium model 
is less bound here and we are not perfectly in the minimum of the jellium deformation 
energy surface, which does not coincide with the deformation associated to the case with 
explicit ions. Consequently, the whole spectra are shifted to lower frequencies where the first 
Iph states are now situated. This corroborates again the fact that the residual interaction 
has a small effect. Moreover, we see by comparison with the Iph spectra that the residual 
interaction can act both ways, repulsive or attractive. It is interesting to note that, in spite 
of the strong general downshift as compared to the case with detailed ionic background, the 
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relation amongst the scissors states (ordering, fragmentation) remains almost unchanged. 
And the downshift is simply given by the spectral gap in the ground state. This hints that 
global deformation is the crucial ingredient determining the pattern of the scissors spectra. 

The comparison of the Iph spectra in the left column of Fig. 3 is puzzling. It looks as if 
the better bound ionic configuration has the higher density of states. Part of this is indeed 
the effect that the larger spectral gap compresses the states just above the gap. Another 
part is generated by the fact that the jellium background is reflection symmetric about all 
three major planes while the ionic background induces slight symmetry breaking. There are 
more degeneracy in the jellium spectrum and none at all for the ionic background. 

The right panels of Fig. 3 show the same analysis for the first isomer. This configuration 



is axially symmetric, i.e. 7 = 0°, but spin polarized with net spin 2 [0. Total deformation 
is (3 = 0.58, close to the deformation of the triaxial state. The upper right panel (axial 
jellium model) shows nicely the exact degeneracy of x with y modes and the absence of 
the z mode. The spectrum is again fragmented in accordance with the high density of Iph 
states. The lower right panel shows results with ionic structure. The z mode still nearly 
disappears. The degeneracy of x and y modes is somewhat broken because a detailed ionic 
structure breaks strict axial symmetry. This cluster is axial only in average, i.e. in the 
sense that 7 = 0°. The Iph spectra are about at the same position in both cases here and 
similarly the fully coupled spectra. The responses are compatible with the assumption of 
small residual interaction. Moreover, we see that the detailed ionic structure does not help 
the z mode to show up (which would, in principle, be possible). Again, we conclude that 
global deformation determines the overall pattern of scissors modes. 

The inserts between upper and lower panels in Fig. 3 show the collective estimates (^ 
and (^ for the scissors frequencies, obtained in a jellium model with given deformation 
parameters. Although the estimates belong to the jellium case, they can be connected with 
lower panels as well, as a collective view where only the global deformation counts. The 
estimates provide rough positions of the scissors modes but they miss any details. The 
actual spectra show that the details, driven by interference with Iph structure, are crucial 
in these very soft clusters with a high level density. 
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B. Naie: crosstalk 

The orbital Ml mode and spin excitations belong to the same channel (see more on spin 
modes in clusters in Refs. ^^-|27]). The orbital strength ~ B(M1) oc N^^^ and power 



spectra ~ B(MlY oc N^^^ evidently overrules the spin strength (oc N^) for large clusters. 
But they can compete in the small clusters considered here. There is an intriguing question 
how strong is a crosstalk between the modes. This point is better discussed in terms of 
power spectra, although we are in the linear regime. Mind that we are now plotting in 
logarithmic scale which is more appropriate for power spectra with their larger variations. 

The various power spectra emerging from excitation of the scissors mode (i.e. with a 
rotational distortion) are shown for Naie in the upper block of Fig. 4 The block exhibits 
angular, spin-dipole and dipole responses. Both unpolarized ground state and spin-polarized 
first isomer of the cluster are considered. This case is different from Nai2 in that both 
configurations are triaxial. The spin-saturated ground state and spin-polarized isomer have 
deformations j3 = 0.35 with 7 = 38" and (3 = 0.32 with 7 = 15°, respectively. 

The uppermost panels show the angular momentum spectra associated with Ml strength. 
The patterns are similar to the previous case of Nai2 but differ in detail due to different 
Iph distributions. The lowest panels of the upper blocks show the power spectrum in the 
dipole channel driven by the crosstalk with the angular momentum mode. The coupling is 
obviously very weak. In fact, there would be no crosstalk at all in a reflection symmetric 
jellium model. The effect here is triggered by the slight asymmetries in the detailed ionic 
configuration. 

The middle panels show the crosstalk with the spin-dipole mode. There is no crosstalk 
at all in the case of the spin saturated ground state. But a strong crosstalk appears for 
the spin-polarized isomer. In that case, the scissors and spin motions are fully coupled and 
have to be considered as one hybrid mode. As complementing information, we show in the 
lower block of Fig. 4 the spin-dipole strength after spin-dipole excitation. For the isomer 
(right column), one recognizes the same spin-dipole spectrum as in middle panel of the 
upper block. The strength distribution is, of course, different because the excitation channel 
was different. For the spin saturated ground state (left column), there is now a clear spin- 
dipole signal. This mode is clearly separated from the scissors, but it resides precisely in the 
same low-energy range as the scissors mode. An experimental discrimination would require 
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to separate the spin and orbital currents. This question has been extensively discussed in 
nuclear physics where the form factor for inelastic electron scattering may give access to the 
separate information (see, e.g., Ref. The problem is much worse for clusters due to very 
low energies involved. The easiest way is to exploit the scaling with system size and to go for 
large clusters where spin-dipole excitations become relatively unimportant (see discussion 
in Ref. §). 

It is worth to emphasize that the cross-talks discussed above connect the modes of the 
opposite space parity. This becomes possible due to mixing configurations of a different 
parity in the detailed ionic background. 



C. Nag^: pure ionic effects 

As we have seen in the examples above, the scissors mode exists if the electron cloud can 
accomplish rotational-like oscillations against a massively deformed background. The mode 
does not work at all in the case of the perfect spherical symmetry. But clusters do have 
ionic structure which breaks spherical symmetry even if the global quadrupole deformation 
vanishes. It is then conceivable that the electron cloud performs rotational-like vibrations 



against the ionic background. Such transversal optical modes exist in the solid state p8 
Their frequency is proportional to momentum as in the acoustic branch. One may expect 
similar low-energy modes in a finite system with zero global deformation as well. 

To check this, we consider the example of Nag"*". Its configuration consists in two rings 
of four ions toped by one single ion. The configuration has (3 = 0.08 and is thus fairly close 
to global spherical symmetry. It is to be noted that the top ion breaks reflection symmetry. 
Figure 5 shows the angular momentum strength in this cluster. There is some strength 
at 1.5-2 eV and another branch at 3.5-4.5 eV. The latter has position and strength much 
similar to the high-frequency branches in the previous examples. (The slight blue-shift is 
caused by positive net charge as compared to the neutral clusters considered before.) The 
branch is again related to the quadrupole plasmon and is predominantly composed of AA/'= 
2 transitions. However, the low-frequency branch has higher frequency and much weaker 
strength than in the previous examples. In fact it is related to AM = 1 transitions. The 
larger frequency is explained by the fact that we deal here with a small magic cluster with 
a big HOMO-LUMO gap. The much lower strength is no surprise. Note that AA/'= 1 
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transitions change space parity and would not be accessed by Ml motion in a jellium model. 
But the detailed ionic background mixes configurations of different parity and so even-parity 
Li operators can generate, though weak, the transitions of AjV —1 kind. This explains why, 
unlike the previous examples, the first branch is weaker than the second one. 

Such motion may be related to the transversal optical mode. But it is not possible to 
track that down unambiguously because in this small cluster the low-frequency Ml mode 
is predominantly of Iph nature. Larger magic clusters are required for such an analysis. 
Altogether, we have seen magnetic Ml strength in a "spherical" cluster. This strength has 
different sources than the scissors mode. It takes place even at zero global deformation and 
has very fittle strength as compared to the genuine scissors modes. 

V. CONCLUSIONS 

We have presented fully microscopic TDLDA calculations for the scissors mode in small 
Na clusters taking into account detailed ionic structure, triaxiality and spin-polarization. 
The results were analyzed in terms of a simple collective (^macroscopic) picture and of more 
microscopic shell-model view of the mode. To single out the effects of ionic background, a 
comparison to jellium cases was done. Clusters with different global deformation (spherical, 
axial and triaxial) were considered. 

The test case of Nan'*' with jellium background displays straightforward and unfrag- 
mented scissors peak. This case was used to study the various contributions to the residual 
interaction defining the actual peak position. It revealed that Coulomb and exchange- 
correlation terms act in opposite directions, resulting finally only in a small collective energy 
shift. The actual sign depends on details, as e.g. interplay with particle-hole structures. The 
calculations for Nai2 show that jellium results demonstrate good agreement with macroscopic 
estimates for the low-energy scissors branch, though the estimates miss such important as- 
pect as fragmentation of the strength due to particle-hole states. There is also a high-energy 
scissors branch which is explained as a result of the transitions through two quantum shells. 

Triaxiality produces a richer picture. It gives rise to an additional L^-scissors branch and 
to considerable increase of the fragmentation. Besides that, triaxiality combined with the 
delicate balance between Coulomb and exchange-correlation contributions can lead to both 
red and blue shifts in one and the same spectrum. 
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Detailed ionic structure leads to the blue-shift of the scissors mode as compared to jellium 
results. It also considerably enhances the fragmentation. The ionic background causes small 
parity-mixing which results in a crosstalk, although weak, between positive-parity scissors 
and negative-parity dipole excitations. This crosstalk takes place for both spin-saturated 
and spin-polarized states. Spin-polarized states demonstrate in addition a strong crosstalk 
of the scissors and spin-dipole modes, the effect which is absent in spin-saturated states. 
Both modes happen to reside in the same corner of low-energy states which makes their 
experimental discrimination extremely tough. 

Comparison of results from jellium with those using detailed ionic background shows that 
the scissors modes arc mainly determined by the global deformation. The ionic background 
as such destroys locally the spherical symmetry and gives rise to some Ml orbital strength 
(as the transversal optical modes in solids). This becomes apparent by the existence of Lj 
response even in clusters with zero global deformation, like Na^ . The strength of such a 
mode is an order of magnitude smaller than the collective scissors strength which lives up 
in truly deformed systems. 

The scissors mode is a general dynamical effect peculiar to different Fermi (and even 
Bose) systems with a disturbed symmetry. In the particular case of atomic clusters, the 
mode demonstrates very interesting features deserving further study. The main problem, 
not tackled in this paper, is its unique experimental identification. Such questions will be 
considered in a forthcoming publication. 
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Figure captions: 

Figure 1: Macroscopic view of scisors mode 0: a) rigid rotation, b) rotation within a 
rigid surface. 

Figure 2: Scissors photo-absorption cross section (upper panel) and radial profile of 
residual interaction (lower panel) for Nan^. Four cases are compared: no residual interaction 
at all (denoted ph, dashed line in upper panel), only Coulomb residual interaction (dotted 
lines), only exchange-correlation part of residual interaction (dashed line in lower panel), 
full residual interaction (full line). No scissors strength could be shown for the exchange- 
correlation part alone because in this case the residual interaction is so strongly attractive 
that an imaginary frequency emerges. 

Figure 3: Strength distribution of angular momentum modes in Nai2. The quadrupole 
deformation of the jellium background is chosen to be the same as in the states with full 
ionic structure. The modes Lx-, Ly, and Lz are distinguished by the lines as indicated in 
the left upper panel. Inserts under the panels show the Iph spectra (bars) up to ~ 2 eV. 
The boxes in the left and right inserts between upper and lower panels show the collective 
estimates (|^) and (^, respectively, for the scissors frequencies in the jellium model. 

Figure 4: The upper block shows various power distributions in Naig after angular ex- 
citation: angular (upper), spin-dipole (middle), dipole (lower). Unpolarized triaxial ground 
state (left) and spin-polarized axial isomer (right) are considered. The lowest block shows 
spin-dipole power following spin-dipole excitation. In all the calculations, full ionic struc- 
ture is taken into account. Every case has three modes in each spatial direction. They are 
distinguished by lines as indicated. 

Figure 5: Strength distribution of angular momentum modes in Nag^, described with 
detailed ionic structure. The Iph spectra at low energies are indicated in the insert below. 
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